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Abstract 

We extend our analysis of the effects of the interplay of quantum phases 
and nonlinearity to address saturation effects in small quantum systems. We 
find that initial phases dramatically control the dependence of self-trapping 
on initial asymmetry of quasiparticle population and can compete or act with 
nonlinearity as well as saturation effects. We find that there is a minimum 
finite saturation value in order to obtain self-trapping that crucially depends 
on the initial quasiparticle phases and present a detailed phase-diagram in 
terms of the control parameters of the system: nonlinearity and saturation. 
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I. INTRODUCTION 



We analyze here some features of a quasiparticle interacting with lattice vibrations that 
are strongly anharmonic in nature. An example of such an oscillator is to be found in 
liquid crystals where the vibration is rotational 0. Work done on these lines by Kenkre 
and collaborators on rotational polarons has shown the existence of fascinating phenomena 
like saturation of selftrapping on increasing nonlinearity beyond a characteristic value . 
Further work on a different kind of oscillator, one which is harmonic for small displacements 
from equilibrium but subjected to a restoring force, which is logarithmic for large displace- 
ments also showed similar behaviour ||. In a study of the more commonly encountered 
harmonic systems, it was shown || that the initial quantum phases can profoundly control 
the process of self-trapping in nonlinear quantum dimers where the interacting vibrations 



were assumed to be harmonic (see also The purpose of this note is to report some 

interesting behaviour that occurs when the initial quantum phases can interact with the 
nonlinearity in the system as well as the saturation effects of anharmonic vibrations. We 
point out here that throughout this paper, we shall assume the vibrations to be classical. It 
is important to note that, although the regime of validity of the DNLSE as a consequence of 
microscopic dynamics is limited [7|-[|, the problem of the interplay between qunatum phases 
and nonlinearity studied in this paper is of general interest in quantum nonlinear equations 
of evolution. The results of the present paper have to be interpreted in this context. We 
emphasize that, throughout this paper, as in earlier work [p|,|To|-[l^| , by the phrase 'quantum 
phases', we refer only to the quantum mechanical phases associated with the quasiparticle. 
Clearly, after the semiclassical approximation has been made, there is no meaning to the 
association of a quantum phase with the vibrational variables. 

We begin here with the following coupled equations for the generalized displacement 
x m of the oscillator at site m, and for the quasiparticle amplitude c m at site m, 

dt ^ / V mn c n -\- .E(x m )c m , (1.1) 



2 



dt 2 



+ u 2 f(x m ) + RE'(x m )\c m \ 2 = 0, 



(1.2) 



where V mn is the intersite matrix element between sites m and n, R is a proportionality 
constant related to the geometry of the system, E( generally nonlinear function 

of the coordinate x m describing the quasiparticle-vibration interaction, and u> 2 f(x m ) is a 
general nonlinear function giving the restoring force on the vibrational coordinate. 

For the rotational polaron, as noted elsewhere 0-fj], the restoring force and the interac- 
tion are described as 

E 



Following the analysis of Kenkre et al. PH, one obtains the amplitude equation 
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dCr, 



X\Cn 



dt V V 1 + (X/A) 2 |c m | 4 

in general, and for the two-site system (dimer) in particular, 

dg{p) 
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p = —2vq , q = 2vp-\ ; — r , r 

dp 



dp 



where 



\ 



g(p) = A 

with the definitions 
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;i.6) 



p = |ci| 2 - |c 2 | 2 ; q = i(c*C2 - cic£) ; r = c*c 2 + C1C2. 



;i.t) 



and where po ; r o refer to the initial values of p,r respectively. Here the nonlinearity \ = 
EqR/uo 2 and the saturation energy A = Eq/A. 

As has been shown in Refs. 0-§J, one obtains from Eq. ( |1.6|) an equation of motion for 
a fictitious conservative oscillator whose displacement is p, and that moves in a potential 
given by 



U(p) = 4VY + g 2 (p) - 2g(p)[g( Po ) - 2Vr ] 



;i.8) 
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II. ROTATIONAL POLARON 



We start this section with the elementary observation that as long as g(p) in Eq. ( |1.8|) 
is an even function of p, the potential will also be one. Furthermore, if g(0) = 0, U(0) = 0. 
This implies that the critical condition to determine the transition between free evolution to 
self-trapped behaviour can be found by examining the condition U(po) < 0. This condition 
combined with Eq. ( |1.8| ) gives 

\g(p Q ) - 2Vr \ >2V. (2.1) 

It is therefore instructive to view the condition expressed in Eq. (|2.1| ) in terms of a phase- 
diagram with the saturation ratio A/V and the critical nonlinearity ratio Xc/V as the 
parameters for different values of po, r marking regions of free and self-trapped evolution. 
We first examine the case for real initial conditions, i.e., q(0) = q = 0. In the next section, 
we shall generalize to complex initial conditions. 

Figure 1 shows the phase diagram for the case when ro < 0. It is to be interpreted in 
the following manner. For a given value of the saturation ratio A/V, as one increases the 
value of nonlinearity ratio Xc/V, one traverses regions of free-evolution. As soon as one 
encounters the boundary (for the relevant value of po) , one enters the region of self-trapped 
evolution. This is the region where Eq. ( |2.1| ) is satisfied. Note that upon increasing the value 
of nonlinearity beyond a certain value one crosses this region and passes into free evolution 
again. This is as expected and is in accordance with the observation in Refs. fUJi]] that 
self-trapping can be destroyed by increasing the nonlinearity. Note also that for the initial 
phase chosen so that r < 0, for a given value of saturation, it requires less nonlinearity to 
enter the self-trapped phase if the amount of initial asymmetry (measured by po) is less. 
This is in consonance with the observation of Ref. ||. It is interesting that this interplay 
between choice of initial phase and initial asymmetry remains qualitatively unaffected even 
in the presence of saturation effects. It is important, also, to note that one gets out of this 
self-trapped phase for correspondingly lower values of nonlinearity. It is interesting to see 
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that as one moves to lower values of po, the system tolerates lesser values of saturation to 
enter the self-trapped phase, the most stringent limit being posed at the value po — 1. It is 
straightforward to estimate the least saturation ratio required (and the corresponding value 



of critical value of nonlinearity ratio): A c /V = ^27/4, Xc/V = y27/2. It occurs for po — > 
and is in agreement with what is obtained in Ref. [||]. There are several features of interest 
in the phase diagram. For fixed p$, the phase boundary exhibits two asymptotes. The first 
of these occurs for finite Xc/V and A/V — > oo. This is the standard DNLSE limit and it 



corresponds to x C A- The value of this asymptote is given by Xc/V = 4(1 — J 1 — pi) /pi 
and is the result of Ref. || . The second asymptote is actually more interesting since it marks 
the boundary marking the destruction of self-trapping. The behaviour along this boundary 
is given by the relationship 

r Al 3/2 



Xc 

V 



V 



a(po), (2.2) 



where a(p ) = p Q y/2/[(l — pfi)(l — yl — pl)} 1 ^ 2 - Equation ( |2.2j ) shows that for small and 
intermediate values of p$, the dependence of Xc/V on A/V is of power-law type with an 
exponent of 1.5. The case of po — 1 has to be dealt with carefully, however, and an analysis 
shows that when p Q — > 1, the equation for the boundary becomes Xc/V = (A/V) 2 . The 
exponent thus changes from 1.5 to 2 and this explains the clear gap in the upper phase 
boundary between small and large values of Pq. 

Let us now examine the case r > 0. Figure 2 is the r > counterpart of Fig. 1. One 
immediately notices the profound difference between the two cases r < and r > 0. For 
instance, one has the largest available area of parameter space for the self-trapped phase 
when po — 1 and this area progressively decreases when p decreases from this value. In 
contrast to the case where r < 0, the phase boundaries for different values of po lie embedded 
within one another in a non-overlapping fashion. It is instructive and straightforward to 
calculate asymptotic values of the various features of the phase diagram especially in the 
region of small po- In this limit, the turning point of the phase boundary (marking the 
mininum amount of staturation allowed for self-trapping) scales with the initial population 



difference as A c /V ~ 6y/3/pQ and the corresponding value of critical nonlinearity is Xc/V — 
\[2A C /V . Again, in this limit, there are two asymptotes for the phase-boundary. The first 
occurs when one is in the standard DNLSE regime and \ <C A. Here, one recovers the 
standard result that the critical value of nonlinearity needed to trap, Xc/V = 8/p^ and 
A/V — > oo. The second asymptote marks the boundary where one gets out of the self- 
trapped phase to the freely evolving phase. This is in the strongly anharmonic phonon 
regime. The relationship between the parameters at this boundary is markedly different 
from the asymptote discussed above and the following equation for the boundary holds: 
Xc/V ~ po(A/V) 3 / 2 . It is important to note that this power-law dependence of Xc/V on 
A/V is the similar to the one found for the case r < 0. However, the crucial point is that 
when r < 0, the dependence of Xc/V and A/V on p Q is quite weak, particularly for p far 
from 1. On the other hand, for ro > 0, both Xc/V and A/V diverge as 1/Pq. This explains 
why, whereas for r < 0, on the upper phase boundary, the curves seem to merge for different 
values of po, the corresponding curves are well-separated for r > 0. There is thus a clear 
quantitative difference between the two quantum phases r < and r > 0. 



III. COMPLEX INITIAL AMPLITUDES 

Thus far, we have only examined the case where the initial conditions for the quasiparticle 
are real, i.e., go = 0. However, it is important to understand the effect of arbitrary initial 
quantum phases. To this end, we write the basic condition for examining the self-trapped- 
free phase boundary (|2.1|) as 

\g(p ) ± 2VyJl-pl-q$\ > 2V (3.1) 

Using the expression ( |1.6| ) in Eq. ( |3.1| ) above, we plot in Fig. 3, the parameter space indexed 
by the saturation ratio A/V and the nonlinearity Xc/V separating the self-trapped regions 
from the free region. For both Fig. 3(a) and Fig. 3(b), r < 0. In Fig. 3(a), po = 0.1 and in 
Fig. 3(b), po = 0.8. The primary observation is that relaxing the constraint of reality of the 
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initial quantum amplitudes results in a loss of parameter space available for self-trapping. 
For small values of po (illustrated by Fig. 3(a)), the phase boundary is much more sensitive 
to changes in go- I n fact, for large and intermediate values of go and for small values of 
Po, the diagram looks very much like the one obtained in Fig. 2 for tq > 0. This can be 
understood very easily by analyzing the expression in conjunction with the equality 
condition of Eq. fl3.1|) that results in 

A = 2(l-|l-pg-gj] 1/2 ) (3 2) 

V ^+l)»+(^) J + ^ta>-l)'+(f) > -2^ + (f) J ' 
For small values of po and nonzero values of go, the right-hand side of Eq. ( |3.2j ) tends to 
diverge as I /pi just as would happen for r > 0, g = because the numerator does not 
vanish in the limit po — > 0. Furthermore, the values for the turning point marking the end 
of the self-trapped region also have the same qualitative behaviour as in the r > 0, go = 
case. Specifically, when g ^ 0, r < 0, the values of the parameters at the turning point are: 
A/V = 3y/3(l - yjl - qD/pl, Xc/V = V2A/V. Thus we see that the values of the turning 
point diverge as 1/Pq, similar to the case when r > 0, g = 0. 

It is possible to understand the qualitative similiarty between the r > 0, g = and 
ro < 0, go 7^ cases. In doing so, we can understand a general similarity between the results 
here and those of Ref. ||. The point is that when r < 0, go ^ 0, in order to trap the 
system (either the DNLSE one or its generalized cousin), one has to bring the system to 
a state that is further away in character from the stationary state of the system, marked 
by r < 0, g = 0, and thus closer to one denoted by r > 0, g = 0. Thus both the cases 
r < 0, go 7^ and r > 0, go = are similar. What is interesting is that this similarity is 
very deep and persists right through into the strongly anharmonic regime. 



IV. LOGARITHMIC OSCILLATOR 



Along lines identical to those discussed above, we have analyzed the free-selftrapped 
phase boundary for the system of a quasiparticle in strong interaction with a logarithmic 
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oscillator discussed earlier in Ref . |5j . The counterpart of Eq. ( |1.6|) is given by 



= ^t e ~x/2xo S inh 2 




(4.1) 



where xo is a saturation parameter. The DNLSE limit is recovered in the limit x/Xo 1. 
The phase-boundary in the parameter space Xo/V,Xc/V marking the transition from free 
evolution to self-trapped behaviour is given then by the expression 



Since the phase-diagrams obtained for this system are qualitatively similar to the ones 
obtained for the rotational polaron, we do not discuss them in detail here. The logarithmic 
oscillator, besides possessing qualitatively features (similar to the rotational polaron) like 
saturation of nonlinearity also exhibits the same asymptotic forms for the phase boundaries 
and turning points. For example, the least saturation value required to support self-trapping 
occurs for r < 0,po — > for Xo/V = e/(2v^2),Xc = 4xo- The asymptotes corresponding 
to Xc/V,Xo/V ~ > °°? are a bit different for p <C 1. In contrast to the power-law type 
behaviour (with exponent 1.5) exhibited for the rotational polaron case, the dependence of 
critical nonlinearity Xc on saturation parameter xo is a bit more complicated and numerical 
evidence seems to indicate the exponent is only slightly greater than 1. However, for po — 1, 
we obtain a result similar to that obtained for the rotational polaron, viz., that Xc/V ~ 
(Xo/V) 2 . For the case, ro > 0, the available parameter space for self-trapping again shrinks 
like the counterpart for the rotational polaron. It is particularly simple to calculate the 
asymptotic properties and turning point of the phase boundary (marking the minimum 
saturation required for self-trapping) in the limit of small values of po- In this regime, 
the turning point is marked by Xo/V — 8 e /Po'Xc = Xo- The lower asymptote (in the 
DNLSE regime) is recovered: Xc/V = 8/pl, Xo/V — > oo whereas the upper asymptote in the 
anharmonic regime is given by the implicit relation: 




(4.2) 



4x c 



(4.3) 



Xo = 



\°g{Xc/V)pl 



,2 ' 



Xc 



OO. 



S 



V. CONCLUSION 



We have extended our analysis || of the effect of quantum-mechanical quasiparticle 
phases on self-trapping by treating saturating dimers . We find that initial quasiparticle 
phases can compete or act with nonlinearity as well as saturation effects. For the case tq < 
(in contrast to r > 0), the region of available parameter space for self-trapping increases 
with decreasing asymmetry of initial population, po- in particular, as reported in Ref. [|J, 
there exists a minimum value of saturation ratio A/V below which no self-trapping can 
occur for any value of nonlinearity. What is interesting, however, is that this critical value 
is the least when tq < 0,po — > and increases for all other regimes. The diametrically 
opposite case of r > 0, go = ®,Po — ► yields a starkly different condition, viz., that this 
critical value of saturation diverges as 1/pq. A generalization to the case of complex initial 
amplitudes gives results similar to the case when r > 0, q = 0. We have also analyzed 
the logarithmic oscillator and obtained qualitatively similar features as those shown by the 
rotational polaron. We are thus led to conjecture that these features are perhaps universal 
to dimers interacting with oscillator potentials that have saturating nonlinearity. 

We gratefully acknowledge V. M. Kenkre for extremely useful discussions. 
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FIGURES 

FIG. 1. The critical value of nonlinearity ratio Xc/V, i.e., the value required for self-trapping 
is plotted logarithmically against the saturation ratio A/V for various values of initial population 
difference po- Here, ro < 0, qo = 0. 

FIG. 2. Xc/V is plotted logarithmically against A/V for different values of po- In this figure, 
r > 0, q = 0. 

FIG. 3. Xc/y is plotted logarithmically against A/V for non-zero values of q^. In (a), p = 0.1 
and in (b), po = 0.8. For both (a) and (b), ro < 0. 
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